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ABSTRACT 

Astrometric detection of a stellar wobble on the plane of the sky will provide us a next break- 
through in searching extrasolar planets. The Space Interferometry Mission (SIM) is expected 
to achieve a high-precision astrometry as accurate as 1 /xas, which is precise enough to discover 
a new-born Jupiter mass planet around a pre-main-sequence (PMS) star in the Taurus- Auriga 
star forming region. PMS stars, however, have circum-stellar disks that may be obstacles to the 
precise measurement of the stellar position. We present results on disk influences to the stellar 
wobble. The density waves excited by a planet move both of the disk's mass center and the 
photo-center. The motion of the disk mass center induces an additional wobble of the stellar 
position, and the motion of the disk photo-center causes a contamination in the measurement 
of the stellar position. We show that the additional stellar motion dynamically caused by the 
disk's gravity is always negligible, but that the contamination of the disk light can interfere with 
the precise measurement of the stellar position, if the planet's mass is smaller than ~ 10 Jupiter 
mass. The motion of the disk photo-center is sensitive to a slight change in the wave pattern 
and the disk properties. Measurements by interferometers are generally insensitive to extended 
sources such as disks. Because of this property SIM will not suffer significant contaminations of 
the disk light, even if the planet's mass is as small as 1 Jupiter mass. 



Subject headings: accretion, 
system: formation 

Introduction 



accretion disks — astrometry — planetary systems: formation — solar 



Since the first discovery by Mayor & Queloz 
(1995), more than 120 extrasolar planets around 
main-sequence stars have been counted in the list 
(Schneider 2004 2 ). Most of them were discovered 
by radial velocity measurement of stars (see, e.g., 
Marcy & Butler 1998; Perryman 2000). One of 
the planets that were already found by the radial 
velocity measurement has been confirmed to have 
transits on the star (Charbonneau et al. 2000), 
and a few planets were originally discovered from 
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their transits (Konacki et al. 2003a, b) . A gravita- 
tional microlensing event was reported that it had 
been caused by a planet orbiting around binary 
stars (Bennett et al. 1999). 

It is considered that astrometric detection of 
wobbling stars will be a strong method for dis- 
covering planets in the near future. A star orbit- 
ing the center of mass of the star-planet system 
draws an ellipse on the plane of the sky (after the 
parallax and the stellar proper motion are sub- 
tracted). The semi-major axis (in arcsec) of the 
ellipse caused by a planet having a circular orbit 
is 



M p r p /AU 
M + M p d/pc 



(1) 



where M p and M are the masses of the planet and 
the star, respectively, r p is the distance between 
the planet and the star, and d is the distance of 
the star from the observer. When a solar mass 
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star with a Jupiter mass planet orbiting at 5 AU 
is observed from 10 pc, the amplitude of the stellar 
wobble is about 500 ^as. 

Measurement of stellar positions with a rela- 
tive accuracy of sub-milliarcsec is a challenge for 
ground-based optical observations. Recent im- 
provement of interferometers has made it possi- 
ble, and planet survey projects are ongoing at the 
Keck Interferometer and the VLTI interferometer 
(for reviews, Perryman 2000; Monnier 2003). 

The Space Interferometry Mission (SIM), which 
is scheduled for launch in 2009, is expected to 
make a significant progress of high-precision as- 
trometry (Danner & Unwin 1999). The SIM is ex- 
pected to achieve ~ 1 /zas relative accuracy with 
its narrow-angle astrometric mode. After the 5 
yr mission duration, it can detect a Jupiter mass 
planet at 300 pc, or a Neptune mass planet at 
20 pc, if it orbits at 1-5 AU around a 1 M Q star 
(Sozzetti et al. 2002). This high accuracy provides 
us an opportunity of searching new-born planets 
in star forming regions. For example, the Taurus- 
Auriga molecular cloud, which is at 140 pc and has 
more than 200 T Tauri stars (Kenyon & Hartmann 
1995), will be a good site for searching a young 
planet. Discovery of a young planet around a pre- 
main-sequence (PMS) star will tell us when planet 
formation occurs. The standard scenario of planet 
formation assumes that planets form through co- 
agulations of planetcsimals (Hayashi, Nakazawa, 
& Nakagawa 1985; Lissauer 1993). According to 
an early version of this model (Pollack et al. 1996), 
it may take more than 10 7 yr to form a Jupiter 
mass planet at 5 AU. However, recent results on 
the oligarchic growth of planet cores (Kokubo & 
Ida 1998, 2000; Thommcs, Duncan, & Levison 
2003) and considerations of dust opacity in a pro- 
toplanet's envelope (Ikoma, Nakazawa, & Emori 
2000; Podolak 2003) indicate that the formation 
of Jupiter may take much less time. In relatively 
massive disks where the asymptotic core masses 
are large, gas giant planets can also form rapidly 
and readily. But in these disks, gravitational in- 
stability may also lead to the rapid formation of 
gas giant (Boss 1997). Although the formation 
time scale cannot be used to directly distinguish 
the paths of planet formation, dynamical diver- 
sity may provide useful clues. In the gravitational 
instability scenario, planets are expected to form 
with large masses and eccentricities. In contrast, 



core accretion is expected to proceed in a tran- 
quil environment with modest mass and low eccen- 
tricity. In the later scenario, planets may acquire 
their dynamical diversity during the disk evolution 
(Ida & Lin 2004) and depletion (Nagasawa, Lin, & 
Ida 2003). A systematic comparison between dy- 
namical properties of protoplanets around PMS 
stars and those of mature planets around main se- 
quence stars provides a useful diagnostic on both 
the formation and evolution of planets and plane- 
tary systems. For example, Boss (1998) and Rice 
et al. (2003) calculated stellar wobbles induced 
by gravitationally unstable disks. Their models 
showed that the disk fragmentation induces sig- 
nificantly larger or faster wobbles than the disks 
without fragmentation. 

PMS stars, however, generally have circum- 
stellar disks. The disk may be a obstacle to the 
planet search using astrometry. Firstly, the planet 
disturbs the disk as well as the star, and moves the 
position of the disk center off the star. The grav- 
ity of the disturbed disk induces an extra motion 
of the star, in addition to the direct wobbling by 
the planet. Secondly, the disturbance of the disk 
breaks the symmetry of the disk emission with re- 
gard to the star. The spiral waves excited by the 
planet make a wavy shape of the disk surface, lead- 
ing to lighting and shadowing the surface. The 
scattered light image of the disk is disturbed so as 
to have a spiral pattern, and the center of the disk 
light is moved away from the star. When the po- 
sition of the star is measured by a small telescope 
(for example, the telescopes of the SIM are 33 cm 
diameter), the beam size is much larger than the 
angular diameter of the disk, thus the observable 
value is only the central position of the light from 
the star-disk system. If the shift of the photo- 
center of the system is similar or larger than the 
amplitude of the dynamical stellar wobble, it is 
difficult to extract the dynamical motion from the 
apparent motion of the photo-center. 

In order to examine how a circum-stellar disk 
affects the precise measurement of the stellar po- 
sition, we study the stellar wobble induced by a 
disk both dynamically (through its mass distribu- 
tion) and illuminationally (through its scattered 
light distribution). In §2, we calculate linear re- 
sponses of the disk to the planet's gravity to obtain 
a spiral wave pattern. For simplicity, the disk is 
treated as two-dimensional and the planet's orbit 
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is circular. In §3, we examine the dynamical ef- 
fect of the disk. The gravity exerted on the star 
by the disturbed disk is calculated, and the stellar 
wobble induced by the disk is compared to that 
by the planet. We show that the disk's dynamical 
effect can be neglected. In §4, the motion of the 
central position of disk light is calculated. We as- 
sume that the three-dimensional wavy pattern of 
the disk surface relates to the surface density fluc- 
tuation calculated in §2, where the disk is treated 
as two-dimensional. The lighting and shadowing 
the disk surface by starlight are calculated to make 
an image of the scattered light. The shift of the 
photo-center is calculated and compared to the dy- 
namical motion of the star. In §5, we apply our 
result to the SIM. We simulate observations by the 
SIM and calculate the anticipated astrometry shift 
induced by the disk light. 

2. Density Wave Excitation by a Planet 
2.1. Wave Equations 

2.1.1. Unperturbed Disk 

We consider a system composed of a central 
star, a circumstellar disk, and a planet in the disk. 
The planet's orbit is assumed to be circular with 
an angular velocity 



0, 



- ( GM 



1/2 



(2) 



where we ignore contributions from the planet's 
mass M p and the disk mass M4, assuming M p , 
M d < M. 

A planet embedded in a gas disk opens a gap 
around its orbit, provided that the planet's mass is 
larger than the Neptune or Saturn's mass [(1 — 3) x 
1CT 4 M ] (Takeuchi, Miyama, & Lin 1996). The 
gap width is in general a function of the planet's 
mass, the gas sound speed, and the gas viscos- 
ity. For a small mass planet like Neptune, the gap 
width is of order of the disk thickness, while for 
a Jupiter mass planet, the gap spreads from the 
m = 2 inner Lindblad resonance (at 0.63r p ) to the 
m = 2 outer Lindblad resonance (at 1.31r p ; Lin 
& Papaloizou 1986; Takeuchi et al. 1996). In this 
work, we consider an unperturbed surface density 
distribution of a disk having a gap, and for sim- 
plicity, it is written as 



co = /g. 



3&*r 



1 

AU 



(3) 



where tau is the non-dimensional radial coordi- 
nate normalized by 1 AU, and 0% is the surface 
density at 1 AU. We assume the surface density 
decreases with r as a power law with an index 1, 
and near the planet's orbit it drops; / gap is unity 
except for the planet's neighborhood at which it 
drops to zero. We use an approximate expression 
for the gap shape: 



/gap — 



exp(±2x) 
exp(±2x) + 1 ' 



(4) 



where the upper and lower signs are for the disk 
outside the planet's orbit and for the inner disk, 
respectively, and the non-dimensional distance x 
from the gap edge is 



x = 



T ^g, out 
^edge 



for the outer disk and 



Aod 



(5) 



(6) 



for the inner disk. The median locations of the 



gap edges 



1 ^(?,out 



and 



are expediently defined 



as the points where the surface density drops to 
half the value of the no-gap disk. The edge width 
Acdgo is of order of the disk half-thickness h. The 
half- width of the gap is A g — (r S;OUt — r 9i i n )/2. We 
assume that the gap shape is symmetrical with 
regard to the planet's position and that A g = 
7"g,out rp — r p rg^iYi. 

A polytropic equation of state with an index 7 
is assumed for the gas. The two-dimensional pres- 
sure is Po = Kgq. The enthalpy and the sound 
speed are written by 



Vo = 



7 



-Kg, 



7-1 



7' 



c o = (7 - l)?7o , 



(7) 



(8) 



respectively. The gas rotation law is assumed 
as Keplerian: £1 2 = GM/r 3 . Though in actual 
disks the gas pressure makes the rotation slightly 
slower than Keplerian, we neglect this small differ- 
ence. For a given enthalpy distribution, the half- 
thickness of the unperturbed disk is h = v / 2?7o/0. 
The derivation of this equation is discussed in 
§4.1.1 below. 
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2.1.2. Perturbation by a Planet 

We solve the gas motions perturbed by the 
planet's gravity. The gap around the planet's orbit 
divides the disk into two parts. The perturbation 
equations are solved separately for the inner and 
outer disks. We use cylindrical coordinates (r, 6). 
For the inner disk the coordinate origin is set on 
the central star, and the coordinates are not an in- 
crtial frame. On the other hand, for the outer disk 
the origin is on the center of mass of the system. 
We neglect the self-gravity of the disk. 

The planet's perturbation potential on the in- 
ner disk is written as 

GM V 



1 ~ (r-2 + r 2 - 2r p r cos(6» - ilpt)) 1 / 2 
GM p r cos(0 - Q p t) 



(9) 



The second term is the indirect term that repre- 
sents the inertial force. For the outer disk, the star 
also orbits around the origin and causes a pertur- 
bation. The perturbation potential is 

GM 



[ r 2 + r 2 _ 2r a rcos(8 - Vt p t - ir)} 1 / 2 



GM n 



[r 2 b +r 2 -2r b rcos(6-n p t)} 1 / 2 ' 



(10) 



where r a = r p M p /(M + M p ) and r b = r p M/(M + 
M p ) are the distances from the origin to the star 
and to the planet, respectively. 

We solve for Fourier components of the pertur- 
bations. Each perturbation variable, which is with 
subscript "1", is decomposed in the azimuthal di- 
rection as 

oo 

Xi(r,0) = Re[X lim (r)exp[im(0 - fl p t)} , 

m=0 

(11) 

where m is the wavenumber in the ^-direction. In 
the following part of the paper, we solve Fourier 
components Xi, m , but we omit subscript "m" in 
writing whenever this omission does not cause any 
confusion. 

The perturbation equations are as follows (Gol- 
dreich & Trcmaine 1979): 

im(Q — Q p )ui — 2Qvi = — ^~(<Ai + ^i) > (12) 



Iffl 

im(n - n p )v 1 + 2Bu x = — — (0i+r?i), (13) 

im(ft-n p )o-i + --^-(ra ui) + — — -v 1 = , (14) 
r dr r 

where u\ and v\ are respectively the radial and az- 
imuthal components of the velocity perturbation, 
and Oort's B = Q + (r/2)(dQ/dr). The enthalpy 
perturbation relates to (J\ as 



2 <7l 

m = c — • 



(15) 



The Fourier transformed perturbation potential 
!>i for the inner disk is expressed as 



Mr) 



GM V (2 - 5 m0 ) 
2r K 



GM ri 



(16) 

where (3 = r/r p , S mn is the Kronecker's delta, and 
is the Laplace coefficient, 

y„(*>-jjf (i-z: < i7 » 

For the outer disk, 

Mr) = (-r- l GM(2 ,: Smo) b? /2 (M 



2r a 

GM P (2 - 8 m0 ) 



2r b 



b? /2 (f3 b ) , (18) 



where f3 a = r/r a and (i b = r/r b . 

2.2. Numerical Integration 

We expect that the planet excites steady waves 
in the disk, i.e., in the frame rotating with the 
planet, the wave pattern keeps the same shape. In 
such a steady state, energy input to the waves from 
the planet is balanced by energy loss of wave dis- 
sipation caused by shocks or the viscous damping. 
This situation is numerically expressed by allow- 
ing the waves to freely escape from the boundaries 
of the computational domain, i.e., we take the ra- 
diative boundary condition. 

We solve the perturbation equations of the in- 
ner and outer disks separately. In solving the in- 
ner disk, the numerical integration starts from the 
point rin i2 near the planet, at which the gas den- 
sity is sufficiently low. Equations (12)-(14) are 
solved inward to a certain radius r in l , where the 
waves wind tightly enough. For the outer disk, the 
computation is performed for [r outil , r outj2 ]. 
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The computational boundaries of planet's side 
are at the gap edges, r; n ,2 = r g - m + 2A d g o for 
the inner disk, and r ou t,i = r g,out — 2A d gc for 
the outer disk. The densities at these boundaries 
are about 2% of the non-gap values (eq. [4]). We 
adopt the free boundary condition. The Lagrange 
derivative of the enthalpy is zero: Drj/Dt = 0. 
Assuming that the radial derivative of the unper- 
turbed enthalpy is zero, we simply have r]i = 0. 

Equations (12)- (14) are integrated using the 
fourth-order Runge-Kutta method. The Laplace 
coefficients (17) are numerically calculated using 
the method by Korycansky & Pollack (1993). The 
boundary condition at n n A or r out ,2 is radiative. 
At these points, only outgoing wave is allowed to 
exist. In a WKB approximation, the enthalpy per- 
turbation is expressed as r]i = H(r)exp[J ikdr], 
where k = ±\J—D/cq and D = ft 2 - m 2 (fl - Q p ) 2 
(eq. [19] in Goldreich & Tremaine (1979) with 
G = 0, uj = mVLp, and k = Q). The radiative 
boundary condition is, as in Korycansky & Pol- 
lack (1993), 



drn _ (}_dH_ 
dr ~ \ H dr 



ik 



(19) 



where we take a positive sign of the wavenumbcr 
k (the negative sign corresponds to the ingoing 
wave) . The derivative of the amplitude of r/i is cal- 
culated by a similar way deriving equation (18b) 
in Goldreich & Tremaine (1979), and is 



1 dH 
H~d7 



ld_ 

2dr 



In 



r<j{)k 
D 



(20) 



For the m = 1 mode in the inner disk, there is no 
wave-zone (i.e., k is pure imaginary). We use the 
free boundary condition (771 = 0) both at n n ,i and 
n n ,2- 

2.3. Viscous Damping of Waves 

The calculations in the previous subsections do 
not take any wave damping into account. In actual 
disks, waves attenuate through viscous dissipation 
or shock dissipation. We consider here viscous 
wave damping, which works even though the wave 
amplitude is not large enough to cause shocks. 
Waves are excited near the planet and propagate 
away. As they propagate further, their amplitudes 
decline through the viscous dissipation. We write 
the enthalpy perturbation in a viscous disk as 



?7l,vis = /vis»7l 



(21) 



where 771 is the enthalpy solved for the non- viscous 
disk in the previous subsection, and /vis, a factor 
of viscous damping, is unity near the planet and 
goes to zero with the distance from the planet. 
From equations (A10) and (A15) in Takeuchi et 
al. (1996), the damping factor is written as 



/ V i s = exp 



-/IX 



4 r> 2 
3 + m 2 (n-ri p y 



m(Q p — il) 



2c 2 



kdr 



(22) 



where r^. m is the radius of the Lindblad resonance 
of the m-th mode wave, and C and v are the bulk 
and shear viscosities, respectively. For simplicity, 
we use the a prescription for viscosity: C = v = 
dvisC^/^t. The radial wavenumber is k = ^/—D/cq 
for outgoing waves. 



2.4. Wave Properties 

In this subsection, we present properties of the 
density waves excited by a planet. The wave pat- 
terns are different from the ones derived by Ko- 
rycansky & Pollack (1993) because of the assump- 
tions that the disk has a gap and that the waves 
viscously dissipate. 

As a fiducial model, we adopt the following pa- 
rameters: The central star's mass is M = IMq. 
The planet's orbital radius is r p = 5AU. The sur- 
face density profile is 00 = 1000/ gap r^u g cm~ 3 

0.4r p and the 



with a gap of the half width A fl 



edge width A 



edge 



5 x 10- 



The disk inner 



radius is r^i = 0.01r p and the outer radius is 
Tout, 2 = 5r p . The polytrope index is 7 = 1.4. The 
sound speed at 1 AU is c (l AU) = 10 5 cm s^ 1 . 
The viscosity parameter is a v i s = 10~ 2 . The 
planet's mass is arbitrary in linear calculations. 
We take a Jupiter mass planet (M p = 1O~ 3 M ) 
as a typical value in the following discussions. 

Figure 1 shows which modes arc effectively ex- 
cited. The angular momentum flux of waves is a 
convenient value to compare, and its expression is 
written by Goldreich & Tremaine (1979) as 



F A 



irmrao 



n 2 -m 2 (n-n p ) 2 

dm 



- Rc(0i + 771 )Im- 



dr 



Im(j7i)Re— (0i + 771) 
dr 

(23) 



where Re and Im mean the real and imaginary 
parts, respectively. Fluxes Fa are evaluated for 
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waves without viscous damping. For such waves, 
Fa converges to a certain value at locations far 
enough from the planet (see Fig. 7 in Korycansky 
& Pollack 1993). In Figure 1, converged values of 
Fa are plotted. Except for small fluctuations, the 
angular momentum flux is a decreasing function 
of mode m, and waves of to < 3 dominate. High 
to > 3 waves are not excited strongly, because 
their Lindblad resonances (LRs) are buried in the 
gap (Artymowicz & Lubow 1994). The locations 
of LRs are closer to the planet for higher m (for a 
schematic illustration, see Fig.l in Takeuchi et al. 
1996). The resonance positions of the to — 3 wave 
are tl,3 = 0.76r p and 1.21r p , which are inside the 
gap of half width A g — 0.4r p , while to = 2 LRs 
at 0.63r p and at 1.31r p are marginally on the gap 
edges. 

Waves cannot freely travel through the disk. If 
the disk viscosity is large enough, waves attenuate 
before they reach the disk inner or outer bound- 
aries. Figure 2 shows the damping factors, / v i s , 
for various to of waves. Higher to waves dissipate 
more quickly and lower m waves can travel further 
away from the planet. If the disk viscosity is large 
(a v i s = 10~ 2 ), however, even m = 2 wave cannot 
reach the central star nor the outer part further 
than 4r p . An m = 2 wave can reach the central 
star if the viscosity is as small as a v i s = 10 -3 , 
although it cannot travel further than 8r p . Even 
in such a low viscous disk, m > 5 waves dissipate 
before they reach the central star. 

The above results are summarized as follows. 
After the opening of a gap by a planet, excita- 
tion of high to waves are suppressed. Only to < 3 
waves can be excited. These waves may propa- 
gate close to the central star, although they dissi- 
pate before arriving at ~ 4r p in a viscous disk of 
a v - ls = 1CP 2 . The wave propagation changes the 
disk shape inside ~ 4r p . Figure 3 shows the sur- 
face density perturbation. The perturbations of 
to = 1 to 20 are summed up to construct the wave 
pattern. In the inner disk, however, the m = 2 
wave dominates the wave pattern, while in the 
outer disk the m — 1 wave dominates. This is 
because there is no inner LR for m = 1, and be- 
cause m > 3 LRs are buried in the gap. The wave 
amplitude declines under the viscous damping at 
larger distances from the planet. 



3. Dynamical Stellar Wobble 

A planet wobbles its central star. While the 
planet's gravity shifts the stellar position directly, 
if there is a circumstellar disk, the planet affects 
the stellar position through altering the disk den- 
sity distribution as well. Since the m = 1 wave 
pattern does not have any symmetry that con- 
serves the center of mass of the disk, the disk's 
mean position moves. Because the mass center of 
a system composed of a star, a planet, and a disk 
must be conserved, the star moves so as to com- 
pensate the motions of the planet and the disk. 

The disk's gravity to the star is considered 
in this section, although the disk self-gravity is 
ignored in the calculation of wave excitation. 
This treatment is justified if the perturbation of 
the disk's gravitational potential is much smaller 
than the gas enthalpy perturbation. When the 
Toomre's Q value, which is cgfl / (ttGoq) , is much 
larger than 1, this condition is satisfied [see also 
the discussion by Ward (1986)]. In our model, 
Q = 19 at 10 AU. 

In the previous sections, different coordinate 
systems have been used for the inner and outer 
disks. We use hereafter the coordinates (r, 8, z) 
and the corresponding Cartesian coordinates 
(x,y, z), whose origin is set to the star. When 
we need to refer to coordinates with the origin on 
the center of mass, the coordinate variables are 
distinguished by writing (r',6', z). 

The x-position of the center of mass of the inner 
disk is given by 

Xd,in = T7 / / J2 Re ^ m ^ 

Md.in y rin ,i JO ~ 

x exp(?TO(0 — Q, p t))]r 2 cos 6d6dr 
= T7 — / r 2 [Re(cr M )cos(fi p i) 

Md.in J rin l 

+ Im(o-i,i)sin(ft p t)]dr , (24) 

where Md; m is the inner disk mass, and the sub- 
script "to" is written explicitly. Because of the 
symmetry, all the contributions vanish after 6- 
integration except for the m = 1 mode. For the 
outer disk the disk variables are calculated in the 
inertial frame, and then 

7j- / >r out,2 

Xd,out = 77 / r /2 [Re(<7i.i)cos(lV) 

M d , out J rLt i 
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+ Im(<7i ; i) sm(Q p t)]dr' + r a cos(Q p t) ,(25) 

where the second term represents the motion of 
the inertial coordinates (/, 8' , z) to the star. The 
y-positions of the centers of mass are similarly 

fin, 2 



Vd.h 



'[Rc((7 U ) sin(ript) 



±v± d.\n J ri n ,l 

lm(a\ t i) cos(Q p t)]dr , 



(26) 



/" r out,2 

Jr't , 



''[Re(a ltl )sm(n p t) 



M d . out 

- Im(<7i,i) cos(Q p t)]dr' + r a sin(Sl p t) .(27) 

The coordinates origin is set to the star, and the 
planet's position is (x p ,y p ) — (r p cos(Q p t), r p sha.(Q p t)) . 
The center of mass of the entire system, which 
does not move in an inertial coordinate system, is 



XCM 



UCM = 



M p X p 
Ms 

M p y p 



YM p 
M A .\ 



■f M d . in 

iUd,in + 



•Md,out2/d,out 



(28) 



(29) 



M s + M p + M d ,in + M d ,out 

The disk contribution to the stellar wobble is 
negligibly small compared to the planet's direct 
effect. Only the m = 1 perturbation shifts the 
center of mass of the disk, and the amplitude of 
the m = 1 wave is not large enough to make a 
significant movement. There is no m = 1 LR in 
the inner disk, and no significant m = 1 wave 
is excited there. For the outer disk the m = 1 
wave is excited at the outer LR (r^i = 1.59r p ), 
which however is not very close to the planet, and 
only a relatively weak wave is excited. Although 
the wave amplitude is proportional to the planet's 
mass, the amplitude of stellar wobble by the planet 
itself is also proportional to the mass. Thus, the 
relative importance of the disk to the planet does 
not change for different planet's masses. 

Figure 4 shows the motion of the star. The 
motion induced both by the planet and by the 
disk (solid line) does not differ significantly from 
the motion induced only by the planet (dashed 
line). In this figure, we assumed that the disk 
mass is 10 times the typical model. The disk mass 
inside 25AU is 0.16M©, which is probably close to 
the upper limit of real disks. Even such a massive 
disk has a minor contribution to the stellar wobble. 
Because the disk contribution is proportional to 
the disk mass, it is negligible for most of actual 
disks. 



4. Wobble of the Photo-Center 

Although a disk does not perturb the stellar po- 
sition significantly, its presence may obstruct the 
precise measurement of the stellar position. In a 
realistic observation, a position of a star is mea- 
sured using a relatively large beam that covers the 
entire star-disk system. The mean location of the 
light from both the star and the disk is interpreted 
as the position of the former. If the disk is dis- 
torted and its light has any asymmetry, it may 
shift the apparent position of the star. In this sec- 
tion, we calculate the motion of the photo-center. 
The photo-center is the mean location of a star- 
disk system that is weighted by its surface bright- 
ness. In our problem, the central star and the disk 
contribute to the photo-center, while the planet's 
radiation is neglected. A Jupiter mass planet is 
1CT 5 Lq at the age of 10 6 yr and 3 x 1CT 6 L Q 
at 10 7 yr, and a 10 Mj planet is ~ 10~ 3 L & at 
the age of 10 6 - 10 7 yr (Burrows et al. 1997). 
The ratio of the luminosities of the planet and the 
star, Lp/L, is much smaller than their mass ratio, 
Mp/M. When a circum-planetary disk forms, its 
emission may contribute (Lubow, Seibert, & Arty- 
mowicz 1999; Tanigawa, & Watanabe 2002; Bate 
et al. 2003; D'Angclo, Hcnning, & Kley 2003). 
We neglect the circum-planetary disk's contribu- 
tion, assuming that most of the disk material has 
accreted to the planet. 

4.1. Scattered Light form the Disk 

In our calculations of the disk radiation, we take 
the following assumptions for simplicity. We con- 
sider optical light images of the disks. At the op- 
tical wavelengths, the disk radiation is dominated 
by the scattered light of the central star, and the 
disk thermal emission is neglected. In addition, 
the disk is treated to be completely opaque. Inci- 
dent starlight cannot penetrate into the disk sur- 
face, but it is completely scattered (i.e., the albedo 
is unity). Thus, when the disk surface waves, only 
the regions that face directly to the star shine, 
making the disk striped. Figure 5 shows a section 
of a disk in the rz-planc. The surface parts that 
are illuminated by the star are indicated by circles. 
The other parts are assumed to be totally dark so 
they do not scatter any light. These assumptions 
make the model disks have more contrasts than 
actual disks. In real disks, starlight can penetrate 
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somewhat inside the surface, and the stripes do 
not have clear boundaries between the light and 
dark parts, making the contrast lower. In addi- 
tion, the assumption that the incident starlight is 
totally scattered at the disk surface without any 
absorption also works to overestimate the disk's 
contrast. Hence, the motions of the photo-centers 
of model disks is more sensitive to the perturba- 
tions than actual. The calculations below should 
be considered to give an upper limit of the photo- 
center shift. 

4-1.1. Variation in the Disk Thickness 

Fluctuation of the disk surface is calculated 
from the results of the wave pattern in the pre- 
vious sections. We assume that the disk always 
keeps hydrostatic equilibrium in the vertical di- 
rection. For the inner disk, 

GMz 



(r 2 + z 2 ) 3 / 2 



GM p z 



[r 2 +r 2 p - 2rr p cos(6> - Q p t) + z 2 ] 3 / 2 



i* =0 

p dz 



(30) 



where p is the gas density, and p is the three- 
dimensional gas pressure. We assume z <C r. 
In addition, outside the Roche lobe (|r — r p \ >• 
(Mp/My/ 3 r p ), the gravity of the planet (the sec- 
ond term) is negligible compared to the gravity of 
the star (the first term). In our model disks, the 
opening gap clears out the gas around the planet's 
orbit, and there remains little gas inside the Roche 
lobe (we ignore the circum-planetary disk). The 
above equation reduces to 



GMz 1 dp 

r 3 p dz 



= 



(31) 



We use the adiabatic equation of state, p = K^p" 13 . 
Then, we have the density profile in the ^-direction 



P = Pmid 



1 



1/(73-1) 



(32) 



where the subscript "mid" means the value at the 
disk midplane, and the disk half-thickness is 



h(r, 0) 




273 Pmid(r,0) 



1 



(33) 



This equation relates the disk half-thickness to the 
midplane density and the pressure, although in 
the previous sections we calculated only the two- 
dimensional values of the density and the pressure. 
It is shown in the Appendix A that a similar rela- 
tion, 



il(r) y 7 — 1 a(r, 0) u(r) 

(34) 

holds for the two-dimensional values. 

The variables of the outer disk are defined on 
the inertial coordinates (r',6', z). The hydrostatic 
equilibrium is expressed as 

GMz 

(r' 2 + r 2 - 2r'r a cos(0 - n p t - tt) + z 2 ) 3 / 2 
GM p z 



[r' 2 + r\ - 2rr b cos(9 - fi p t) + z 2 ] 3 / 2 
1 dp 



p dz 



= . 



(35) 



For the outer disk, r' 3> r a . We assume 
r' ;» z, then outside the Roche lobe (\r' — rb\ S> 
(Mp/M) 1 / 3 ^), the above equation reduces to 



GMz 1 dp 

r- = 



r /3 



p dz 



(36) 



Then, we have an expression of h similar to that 
for the inner disk, 



(37) 

where n 2 (r') = GM/r' 3 . 

A disk height variation induced by a 10 Jupiter 
mass planet is shown in Figure 5. We assume that 
the disk has an opaque surface at h, and starlight 
is totally scattered at the surface. Only the parts 
that are illuminated directly by the star emit ra- 
diation, while the other parts are totally dark. 

4-1.2. Scattered light flux from the disk 

We neglect the star's diameter. The starlight 
emerges from the origin. If there is no obstacle 
on the ray to a part of the disk surface, that part 
is illuminated by the star. The parts marked by 
circles in Figure 5 are illuminated. The starlight 
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is assumed to be totally scattered at the disk sur- 
face (r, 6, h). Though the actual scattering occurs 
at the photo-surface of constant optical depth, 
which is below the disk surface, we approximate 
the photo-surface by the disk surface. In the Ap- 
pendix B, this approximation is examined. Con- 
sider an illuminated part at (r, 9, h) on the in- 
ner disk. An area rdrdO on the xy-plane (mid- 
plane) corresponds to an area of the disk surface 
dS = rdrdO/cosip, where (p is an angle between 
the normal vector n of the surface and the z-axis. 
The surface dS receives a stellar energy flux 

r cos adS 

where L is the stellar luminosity and a is the 
angle between the normal vector n and the line 
from dS to the star. To derive cos a, we assume 
that the waves wind so tightly that dh/d9 is ne- 
glected comparing to dh/dr, i.e., the toroidal com- 
ponent of the normal vector n is neglected and n 
= (n r ,ng 7 n z ) w (— sintp, 0, cos ip). We then have 

cos a = sin (3 sin <p> — cos (3 cos ip , (39) 

where (3 is the angle between the ray of starlight 
to dS and the z-axis. For simplicity, we assume 
the scattering at the disk surface is isotropic. The 
intensity of the scattered light is 

_ dL d _ £cosa 

~ ndS ~ 4ttV 2 ' 1 j 

An observer looks at the disk with a viewing angle 
v (v = at face-on), using a detector that has a 
solid angle dil from the disk. We define the angle 6 
between the line of sight (pointing to the observer) 
and the normal vector n, which satisfies 

cos S = — sin v sin <p cos 9 + cos v cos <p . (41 ) 

The scattered light flux from the surface element 
dS to the detector is then 



dF d 



I cos SdSdfl 
L 



4 7r 2 r 2 



dh 



-1/2 



dh h 
dr r 



x ( cos v — — sin v cos 6 ) rdrdOdQ. , (42) 



where we used 
r 



sin/3 : 



( r 2 + h 2y/2 



1 , cos (i — 



h _ h 

(r 2 + h 2 ) 1 / 2 ~ 7 
(43) 



simp = 



dh 
dr 



1 + 



dh\ 2 _ (dh 

Tr)\ ' C0S " = [ 1+ U 

(44) 

and h-^r. 

Next, consider an illuminated part at (r',9',h) 
on the outer disk. The disk surface dS' = 
r' dr' d9' / costp' receives a stellar energy flux 



-1/2 



dL d = L 



cos a'dS' 
Airr 2 



(45) 



where r is the distance between the disk surface 
element dS' and the star, and 



cos a' = sin (3 sin <p' — cos (3 cos <p' + O 



sin/3 



(r 2 + h 2 ) 1 / 2 



1 , cos (3 



(- 



h 



(46) 



, dh 



-1/2 



( r 2 + h 2y/2 r 

(47) 



C0S</? = 



, dh 
' + 1 d? 



The intensity of the scattered light from dS' is 

dL d L cos a' 
= ndS' = 4n 2 r 2 ' 

The scattered light flux is 
dF d = I cos S'dS'dSl 



(48) 
is 

(49) 



, dh 
: ~ 1 d? 



-1/2 



dh h 
dr' r 



x ^cos v - sin v cos 0'^ r'dr'd9'dfl ,(50) 

where 

cos 5' = — sinw sin 93' cos9' + cosw cos</?' . (51) 

A simulated image of a face-on disk is shown in 
Figure 6. 

4.2. Motion of the Photo-Center 

4-2.1. Face- on Disks 

If a disk is observed at face-on (v = 0) , then the 
most important mode in moving the photo-center 
is the m = 1 wave, which also shifts the mass cen- 
ter of the disk. Modes of m > 2 can also shift the 
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photo-center slightly, because the disk emission is 
a non-linear function of the perturbations (see eqs. 
[42] and [50] ) and does not have a symmetry after 
summing up m > 2 perturbations. The contribu- 
tions of m > 2 modes are, however, small. 

Figure 7a shows the motion of the photo- 
center of the star-disk system. The disk scat- 
tered light significantly affects the motion of the 
photo-center, which is shown by the solid line. 
The motion of the disk partly cancels the stel- 
lar motion, making the amplitude of the photo- 
center's motion about two-third of that without 
the disk (shown by the dashed line). The disk 
contribution to the motion of the photo-center is 
significant. In this figure the planet's mass is a 
Jupiter mass (Mj w 10 _3 M Q ). For a larger mass 
planet, however, the disk's relative contribution 
becomes smaller. Figure 7 b shows that for a 10 
Mj planet the motion of the photo-center with 
a disk does not significantly differ from the mo- 
tion without the disk. The amplitude of the star's 
dynamical motion is proportional to the planet's 
mass, while the disk contribution to the motion 
of the photo-center is not proportional to it. The 
contrast of the spiral pattern shown in Figure 6 
is nearly saturated even though the planet is as 
small as Jupiter. Thus, a further increase in the 
planet's mass does not cause a significant increase 
in the spiral pattern contrast nor a much larger 
shift of the photo-center by the disk. A more mas- 
sive planet results in a decreasing importance of 
the disk in the photo-center's motion. 

4-2.2. Inclined Disks 

Waves of m = 2 and higher do not move the 
center of mass of a disk, and only slightly shift the 
photo-center of a face-on disk. If a disk is observed 
at an inclined angle, however, m > 2 waves also 
contribute in moving the photo-center. When a 
flared disk is inclined to the observer, the half part 
on the farther side is illuminated more efficiently 
than the other half. This imbalance breaks the 
symmetry of m > 2 waves. 

Figure 8 shows the scattered light image of a 
disk that is 75° inclined to the line of sight. The 
bright part of the disk is the farther side to the ob- 
server, and the closer side is dark. The position of 
the photo-center shifts from the star to the bright 
side. This shift is as large as the disk radius multi- 
plied by the disk-to-star luminosity ratio, rdLd/L, 



which is much larger than the amplitude of the 
stellar wobble. In observations of the stellar posi- 
tion, however, it leads to only an additional offset 
to the mean stellar position, because observational 
beam sizes are usually larger than the disk radius. 
An important observational quantity is the mo- 
tion of the photo-center with time. As the planet 
orbits, the wave pattern rotates. This causes a 
fluctuation of the photo-center, because only a half 
disk is on the bright side. Figure 9a shows the mo- 
tion of the photo-center caused by a 1 Mj planet. 
The disk is 60° inclined to the line of sight. We 
show the relative motions from the average posi- 
tion. The asymmetry of disk light between the 
bright part and the dark part causes a large offset 
of the photo-center position from the stellar posi- 
tion, but this offset is subtracted in plotting Figure 
9. Presence of the disk changes the orbital shape 
of the photo-center in addition to enlarging the 
motion. The amplitude of the disk effect is sim- 
ilar to that in the face-on case. Hence, the disk 
light interferes with precise determination of the 
planet's orbit. For more massive planets, the rela- 
tive importance of a disk becomes smaller. Figure 
96 shows that if the planet's mass is as large as 10 
Mj, the disk effect is not very significant. 

4.3. Various Models 

4-3.1. Low Viscosity 

In our typical models, we used a v - ls = 10 -2 
for the viscosity parameter. Here, we consider 
a model with a lower viscosity a v i s — 10~ 3 . As 
shown in Figure 2, waves have attenuated before 
arriving at r out .2 = 5r p for an a v i s = 10~ 2 disk. In 
a low viscous disk, however, waves can reach as far 
as 10r p , so we use a disk model with r out .2 = 10r p . 
Figure 10 shows the motion of the photo-center 
of a 60° inclined disk. The motion of the photo- 
center has a different orbital shape from that of 
the high viscous disk, though their amplitudes are 
similar. The wave spiral pattern in a lower vis- 
cous disk extends to a larger region of the disk, 
leading to a different shift of the photo-center of 
the disk. This result means that the motion of the 
disk's photo-center is sensitive to the wave pat- 
tern, which changes according to the disk proper- 
ties. Thus, if the disk contamination to the photo- 
center is significant, it is difficult to remove this 
contamination. 
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4-3.2. Narrow Gap 

In our typical models, we assume a gap between 
[0.6r p , lAr p ] (where the density is below half of the 
value without the planet). If the planet's mass is 
larger than the Jupiter mass, it clears the disk 
gas from the m — 2 inner LR at 0.63r p to the 
m = 2 outer LR at 1.31r p (Lin & Papaloizou 
1986; Takeuchi ct al. 1996). Thus, the gap size as- 
sumed here is the typical value, but we also made 
an experiment to see what occurs if the gap is 
narrower. The motion of the photo-center in Fig- 
ure 11 is calculated with an assumed gap between 
[0.8r p , 1.2r p ]. The narrow gap allows wave excita- 
tion of modes as high as m « 6, while in the wide 
gap case only m < 3 waves are excited. It causes a 
stronger disturbance in the disk. The fluctuation 
of the photo-center induced by the disk is so large 
that the precise measurement of the planet's orbit 
becomes difficult. 

5. Astrometry Shift Observed by the SIM 

The astrometric wobble in the position of a 1 
Mq star at a distance of 140 pc due to a Jupiter 
mass planet in a 5 AU orbit is ~ 36 //as (~ 7 //as 
for 1 AU orbit). The intcrferometric observations 
by SIM will have the accuracies to measure such 
small angular displacements. However, SIM delay 
measurements correspond to the position of the 
photo-center. In order to derive the astrometry 
parameters of the target star, such as parallax, 
proper motion, and star's wobble about the cen- 
ter of mass, we have to assume that the measured 
photo-center position shifts represent the true as- 
trometry signal of the target. Often it is assumed 
that the target has only non-luminous stellar or 
sub stellar companions. If the target star is not in 
a clean environment, but embedded in a disk, it is 
likely that the position of the photo-center as mea- 
sured by SIM delay is not a true representation of 
the target star position. 

SIM will have the ability to measure the phase 
of the complex visibility of the star with unprece- 
dented accuracy. In the visible wavelengths the 
star is a million time brighter than the planets. 
Therefore, in the idealized case of a single star 
with planets, the brightness distribution can be as- 
sumed to be point source with the photo-center co- 
incident with the star position. The complex visi- 
bility of the brightness distribution has an ampli- 



tude unity with the phase changing as the photo- 
center wobbles around the star-planet center of 
mass. The visibility phase shift A(f> measured by 
the interferometer with a baseline b = 10 m is re- 
lated to the angular separation AO between the 
center of mass and the star position: 

Acb = ^A9 , (52) 
A 

where A is the wavelength. The measured visi- 
bility phase A</> can result from (i) the systemic 
motion of the target star, (ii) the fringe measure- 
ment errors due to the photon noise, and (iii) from 
the brightness distribution around the target star, 
such as the scattered light from the disk. The 
photon and instrumental noises introduce a ran- 
domly varying position shift which can be reduced 
by longer observing time. However, the photo- 
center shift induced due to the brightness struc- 
tures around the star introduces a systematic po- 
sition error, which vary with time only as the ori- 
entation of the baseline changes. This makes in- 
terpretation of the measured astrometry signal in 
terms of wobble induced by planet(s) difficult. It 
can be partly removed by modeling if we know the 
target environment well enough. 

The disk's contamination to the measured 
astrometric shift is especially significant if the 
planet's mass is smaller than about 10 Jupiter 
mass. An interferometer is not too sensitive to 
structures on a much larger spatial scale than the 
fringe size A/6. In the case of an extended source 
with a uniform brightness the interferometer fringe 
visibility decreases with the increasing source size 
a src as sinc(7ra src A/6), where sine x is (sinx)/x. 
It becomes zero when the source size is equal to 
the fringe size. For SIM the fringe size is about 10 
mas at 0.5//m. Although an extended disk with a 
radially smooth symmetric distribution (with size 
3> 10 mas) will be well within the field view of 
SIM, its fringes will be washed out significantly 
and therefore, will not contribute to the fringe 
visibility of the target star. On the other hand, 
any small scale structures (size < 10 mas) present 
in the scattered light within the SIM field of view, 
such as those caused by disk spirals or gaps, can 
produce fringes. Thus, the small scale structures 
in the disk can contribute to the complex fringe 
visibility of the target star, resulting in a phase 
shift that will be interpreted as astrometry signal. 
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In this section, we present the results of the 
simulated astrometry shifts for some of the disk 
light models. In order to simplify the computa- 
tions we assumed 8 spectral channels covering the 
wavelength range 400 to 900 nm. The complex vis- 
ibility phases of the star-disk structure were com- 
puted for each channel for given orientation of the 
baseline. We use SIM parameters b = 10 m and 
mirror size 33 cm. Using equation (52) we esti- 
mate the astrometric position shift for each chan- 
nel. The final values for astrometry shifts were 
obtained by averaging over the channels weighted 
by SIM sensitivity in each channel and by the 
stellar photon flux in each channel. The stel- 
lar magnitudes, spectral type, and estimated SIM 
throughput were used to compute the number of 
photons detected in each spectral channel. For 
SIM throughput we adopted the values obtained 
from the design team for the aperture, transmis- 
sion, and detector efficiencies. We also assumed a 
field stop of 1", but the band width smearing re- 
sulting from combining the astrometry shifts over 
all channels reduces significantly the effective field 
of view to 30 mas. 

For these simulations we assumed the objects 
are at Taurus distance, 140 pc. A Jupiter mass 
planet orbits at 5 AU from a solar-mass star, lead- 
ing to a motion of the star with an amplitude 
36 /xas. The stellar photospheric emission is ap- 
proximated by a blackbody of 4000 K. The disk 
scattered light image is constructed on the nu- 
merical grids of 2001 x 2001. The star is set on 
the central grid at (1001, 1001). Figure 12 shows 
simulated astrometry shifts for the SIM. We show 
only the additional shift caused by a disk, by sub- 
tracting dynamical shifts induced by a planet. A 
disk of our typical model is viewed at face-on. The 
amplitude of the measured variation in the photo- 
center's position is less than 1 /ias, which is much 
smaller than the amplitude of the stellar motion, 
36 ^as. Figure 13 shows the astrometry shifts for 
an inclined disk. The disk is observed at an in- 
clination v = 60°. We show the shifts measured 
with two baselines of the SIM. One is parallel to 
the major axis of the disk image, and the other 
is perpendicular. The amplitudes of the measured 
variations in the photo-center's position are less 
than 1 ^as. Hence, presence of a disk does not 
make a considerable contamination on the mea- 
surement of the stellar position. Figure 14 shows 



a case in which a planet orbits at 10 AU, which 
is twice of the typical models. The disk size is set 
twice, 50 AU. The size of wave pattern also be- 
comes twice, and we expect the amplitude of the 
photo-center's motion is two times larger. The 
predicted astrometry shifts are however compara- 
ble to the case of a planet at 5 AU. The scale 
length of the wave pattern is two times larger, and 
because a larger structure is less sensitive to in- 
terferometer observations, the disk has a smaller 
effect on the astrometry shift. This decrease in 
the sensitivity nearly cancels the increase in the 
photo-center's motion. The dynamical shift by a 
planet is proportional to the planet's orbital ra- 
dius. For a planet at 10 AU, the dynamical shift 
has an amplitude 71 ^as. Hence, the measure- 
ment of a planet with a lager orbital radius suffers 
less contamination of the disk light, although we 
need a longer observing time to detect such a long- 
period planet. 

6. Conclusions 

The influence of a circumstellar disk on the mo- 
tion of the star induced by a planet is analyzed. 
We show the following. 

(1) The m = 1 disturbance in a disk dynam- 
ically moves the stellar position. The disk dy- 
namical effect is, however, much smaller than the 
planet's gravity, and can be neglected. 

(2) Disk emission moves the photo-center of a 
star-disk system. In some cases, the disk consider- 
ably shifts the photo-center and affects the precise 
measurement of the stellar position. A significant 
contamination from the disk occurs if the planet's 
mass is smaller than 10 Jupiter mass, or if the gap 
in the disk is narrow. For more massive planets, 
the disk contamination is less important. 

(3) Actual measurement of the stellar position 
using an interferometer arc insensitive to the ex- 
tended large scale structures in the disk surface 
brightnesses. The SIM measurement, for example, 
does not suffer a significant contamination from 
the scattered light of the disk. 

We consider that the above results provides an 
upper limit of the disk light contamination, be- 
cause our models assume maximum contrast of the 
disk spiral pattern (see §4). In addition, the dust 
residing above the disk scatters the starlight and 
illuminates the disk surface with a large incident 
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angle. This secondary illumination from the upper 
dust also reduces the disk's contrast. Hence, in ac- 
tual disks, the disk light contamination is probably 
smaller. We need a more realistic radiative trans- 
fer calculation in order to take account of such 
effects and to evaluate precise values of the disk 
light contamination. 

In this paper, the planet's orbit is assumed to 
be circular. A planet having an eccentric orbit 
opens a wider gap than a planet of a circular or- 
bit (Artymowicz & Lubow 1994). Because a wider 
gap suppresses the wave excitation, the disk's spi- 
ral pattern is weaker. Hence, for a planet of an 
eccentric orbit, the disk's contamination is less im- 
portant. If the planet's orbit is inclined to the 
disk, however, the planet excites bending waves 
(Shu, Cuzzi, & Lissauer 1983; Artymowicz 1994; 
Ostriker 1994), which may have a large amplitude 
of wavy pattern on the disk surface. In such cases, 
the disk's contamination may be more significant 
than that estimated in this paper, and should be 
examined as a future work. 
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A. Derivation of the Disk Thickness from the Two-dimensional Density and Pressure 

In this appendix, wc derive the disk half-thickness from the two-dimensional values of the density and the 
pressure (eq. [34]). Integrating the gas density p (eq. [32]) in the z-dircction gives the surface density, 

(J = 2 J Q pdz = Pmid hB Q,— + • ( A1 ) 

The beta function B is 



where T is the gamma function. The two-dimensional pressure is similarly 

2' 73-1 



p = 2 J o P dz = PmidhB ( I, _ | ). (A3) 



Using equations (Al) and (A3) with the relation 



B(i^T + l)=^ 1 B( 1 -,^- + l) , (A4) 
\2 73-1 J 373-I V 2 73 - 1 



the disk half-thickness (eq. [33]) is 



,1 / 2 73 p mid _ 1 / 3 73 -I P 

/l ~0V73-lp m id~0V 73-1 a ' (A5) 

From equations (Al), (A3), and h p^^j ^ ; we see the relation between the two- and three-dimensional 
7's as 7 = (373 — l)/(73 + 1). Then, the disk half-thickness is 



which is identical to equation (34). 



B. Location of the Photo-surface 

In §4.1, it is assumed that the starlight scatters at the disk surface (r, 8, h) at which the gas density drops 
to zero. The realistic scattering, however, occurs at the photo-surface of constant optical depth, which is 
below the disk surface. In this appendix, we examine this assumption. The optical depth at (r, 6, z) to the 
starlight is 

-(h—z)p(z)n 



cos a 

1 /h — ?\ T3/(73-l) 
]_ 2 l/(73-D Kpmid/l (^i) (Bl) 

cos a \ h J 

where a is the incident angle of the starlight, k is the opacity, which is ~ 100 for the visible light if the 
dust-to-gas ratio is similar to the solar abundance and the dust is well mixed to the gas (Miyake & Nakagawa 
1993). In the second equality, equation (32) is used. The condition of the photo-surface, r = 1, reduces to 

^ ~ 2-v.s (j^y-- 1 '/" . (B2) 

h \Kp mid hJ 
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We use the values at 10 AU: 73 = 1.4, cos a w (dh/dr — h/r) ~ 0.1, k ~ 100 cm 2 g _1 , /9 m id <~ 10 _11 g cm~ 3 , 
and ft ~ 10 13 cm. Equation (B2) becomes (ft — z)/h ~ 0.02. Thus, the photo-surface almost coincides with 
the disk surface. However, we are considering an era after the formation of gas giant planets. If most of 
the dust has become large bodies, the opacity n is greatly reduced. In addition, the dust sedimentation 
to the midplane also reduces the opacity near the disk surface. These effects may shift the photo-surface 
significantly below the disk surface, and are the subject of future investigations. 
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Fig. 1. — Angular momentum fluxes carried by 
the waves of various modes m. The diamonds (O) 
and crosses (+) show fluxes in the outer disk and 
in the inner disk, respectively. 
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Fig. 2. — Wave damping factor / v ; s . The solid, 
dotted, and dashed lines correspond to m = 1, 
2, and 5 waves, respectively. The damping factor 
(eq. [22]) is calculated for the inner part inside 
the inner Lindblad resonances (ILRs) and for the 
outer part outside the outer Lindblad resonances 
(OLRs). Between the ILR and the OLR, / v j s is 
assumed to be unity. For the m = 1 wave, there 
is no ILR and / v ; s is set to unity in the inner disk 
inside the planet's orbit, (a) a v j s = 1(T 2 . (b) 
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Fig. 3. — The surface density perturbation nor- 
malized by the unperturbed value, ai(r)/<r (r), is 
plotted. A planet is at (5,0), and opens a gap. 



Fig. 5. — Fluctuation of the disk surface. A vari- 
ation in the disk half-thickness in the planet's di- 
rection (9 = Qpt) is shown by the solid line. The 
portions marked by the circles are directly illumi- 
nated by the star. To show the fluctuation clearly, 
we use a model with a 10 Jupiter mass planet 
(M p = 10Mj) at 5 AU. A gap opens around the 
planet's orbit. 
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Fig. 4. — Movement of the stellar position in the 
inertial coordinates. A Jupiter mass planet or- 
bits at 5 AU. The movement caused both by the 
planet and the disk is shown by the solid line. The 
dashed line shows the movement caused only by 
the planet, ignoring the disk's effect. The dotted 
line shows the contribution of the disk (the differ- 
ence between the solid and dashed lines). 



Fig. 6. — Simulated scattered light image of a face- 
on disk. The disk is perturbed by a Jupiter mass 
planet (M p — Mj). Gray scale shows logarithmic 
values of the disk flux. 
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Fig. 8. — Simulated scattered light image of a in- 
clined disk. The disk is perturbed by a Jupiter 
mass planet (M p = Mj) and viewed at v = 75° 
(v = for a face-on disk). The top half of the 
disk is the farther side to the observer, while the 
bottom half is the closer side. Gray scale shows 
logarithmic values of the disk flux. 



Fig. 7. — The shift of the photo-center for a face- 
on disk. The solid line show the motion of the 
photo-center of the star-disk system. The dashed 
line show the motion of the star only by the planet, 
assuming no disk. The dotted line shows the disk's 
contribution (the difference of the above two re- 
sults). The diamonds show the positions of the 
photo-centers when the orbital phase of the planet 
is zero, (a) The photo-center motions with a 
Jupiter mass planet (M p = Mj). (b) The motions 
with a 10 Jupiter mass planet (M p = lQMj). 
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Fig. 9. — The shift of the photo-center. The 
disk is viewed at an inclined angle v = 60°. The 
solid, dashed and dotted lines represent shifts in 
photo-centers (see Fig. 7) (a) The photo-center 
motions with a Jupiter mass planet (M p = Mj). 
(b) The motions with a 10 Jupiter mass planet 
(M p = IQMj). 



Fig. 10. — The shift of the photo-center with a 
low viscous disk a v i s = 10~ 3 . The disk is viewed 
at an inclined angle v = 60°. The solid, dashed 
and dotted lines represent shifts in photo-centers 
(see Fig. 7) (a) The photo-center motions with a 
Jupiter mass planet (M p = Mj). (b) The motions 
with a 10 Jupiter mass planet (M p — 10Mj). 
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Fig. 12. — The predicted astronomy shift by the 
SIM is plotted against the orbital phase of the 
planet. The disk is 140 pc away and observed 
at face-on. A Jupiter mass planet orbits at 5 AU 
from a solar mass star. The planet's orbital phase 
is measured from the baseline position angle of the 
SIM. 
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Fig. 11. — The shift of the photo-center with a 
narrow gap between [0.8r p , 1.2r p ]. The disk is 
viewed at an inclined angle v = 60°. The solid, 
dashed and dotted lines represent shifts in photo- 
centers (see Fig. 7) (a) The photo-center mo- 
tions with a Jupiter mass planet (M p — Mj). 
(b) The motions with a 10 Jupiter mass planet 
{M p = 10Mj). 



Fig. 13. — The predicted astronomy shift by the 
SIM is plotted against the orbital phase of the 
planet. The disk is observed at an inclination 
v = 60°. The solid line shows the shift that is 
measured with a baseline parallel to the minor 
axis of the disk image, while the dashed line cor- 
responds to a baseline parallel to the major axis. 
The planet's orbital phase is measured from the 
minor axis of the disk image. The disk is 140 pc 
away and a Jupiter mass planet orbits at 5 AU 
from a solar mass star. 
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Fig. 14. — Same as Fig. 13, but for a Jupiter mass 
planet at 10 AU from a solar mass star. 
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